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ABSTRACT 

We present spectral analysis of ASCA data of 17 elliptical and 3 lenticular galaxies 
most of which have high L^/L-q- Single-temperature models (MEKAL and Raymond- 
Smith) give unacceptable fits (Xred ^ ^■^) ™ most cases and, in agreement with 
previous studies, give very sub-solar abundances, {Z) = 0.19 ± 0.12 Zq (MEKAL). 
The spectra for approximately half the sample are better fit with a cooling-flow model 
which in three cases gives a substantially better fit. The abundances derived from the 
cooling-flow model are also significantly larger, {Z) = 0.6 ± 0.5 Zq. We empirically 
tested the reliability of the plasma codes in the Fe-L region and found no evidence for 
serious problems with the determined temperatures and abundances. 

Two-temperature models give substantially better fits that are formally acceptable 



red 



1 INTRODUCTION 



The me tallicity of the hot X-ray emitting gas in ellipticals 
is a sensitive diagnostic of their supernova history and thus 
serves as a powerful discriminant for models of the enrich- 
ment of the hot gas (e.g. Loewenstein & Mathews 1991; 
Ciotti et al. 1991). Previous ROSAT and ASCA studies 
have found the metallicities of ellipticals to be very sub- 
solar in contradiction to the standard enrichment models 
( e.g. Sarazin 1997). T his discrepancy has led some authors 
( Arimoto et al. 1997 ) to question the reliability of the stan- 
dard plasma codes, in particular their modeling of the Fe-L 
emission lines. 

However, the abundances derived from the plasma codes 
also depend on the temperature structure of the gas. Be- 
cause the cooling times are short (~ lO'"* yr) in ellipticals 
(e.g. Thomas et al. 1986) the hot gas in elliptical galaxies 
may consist of a continuum of hot and cold phases similar 
to that expected for galaxy clusters (Nulsen 1986; Thomas, 
Fabian, & Nulsen 1987). Fitting single-temperature models 



1.0) in all but a few cases. The highest S/N galaxies (which also have highest 
Lyi/LB) have fitted temperatures < 2 keV for both components consistent with each 
being distinct phases of hot gas. The lowest S/N galaxies (which also have lowest 
L^c/Lb) generally have a hot component with temperature, Th 5 keV, which is 
consistent with emission from discrete sources. (We discuss the origin of these two 
components from analysis of — Lb and L^/Lb-) The abundances of these two- 
temperature models are approximately solar, (Z) — 0.9±0.7Zq (MEKAL), consistent 
with a recent multi-phase model for the evolution of hot gas in ellipticals. Finally, for 
several galaxies we find evidence for absorption in excess of the Galactic value and 
discuss its implications. 

Key words: galaxies: general - galaxies: evolution - X-rays: galaxies. 



to intrinsically two-temperature spectra can result in a sig- 
nifica nt underestimate of the metallicity ( Buote & Canizares 



1991| ). 

Matsumoto et al. (1997) have obtained very sub-solar 
abundances from fitting two-component models to ASCA 
SIS and GIS data of a sample of 12 ellipticals, but they 
restricted the form of the high-temperature component to 
have the same temperature for each galaxy. That is, Mat- 
sumoto et al. combined the data at high energies, _E ~ 4 — 10 
keV, for the 12 galaxies and found that the composite spec- 
trum is not tightly constrained but equally well described 
by either a thermal bremsstrahlung model with T ~ 12 keV 
or a power-law model with photon index 1.8. These mod- 
els, which are thought to describe the emission from discrete 
sources, are then assumed to be similar for every galaxy and 
to apply over the full ASCA energy band, E = 0.4-10 keV. 
However, these models need not be valid representations for 
every galaxy, especially for energies below 4 keV where they 
are only extrapolations. 

Hence, there is a need for a study of a larger sample of 
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Table 1. SIS vs GIS 



Energy 


Area ( 




AE (eV) 


(keV) 


SIS 


GIS 


SIS GIS 


0.75 


60 


6 


70 ~ 200 


1.00 


140 


50 


75 ~ 200 


1.25 


175 


100 


80 225 


1.50 


200 


130 


85 240 


2.00 


150 


160 


90 250 


5.00 


130 


120 


150 400 



Comparison of the effective area and energy resolution of the SIS 
and GIS. These values are approximate as they were obtained 
from visual inspection of Figures 8.4a, 8.4b, and 9.5 of the ASCA 
Technical Description (Appendix E, AN 97-OSS-02). The spectral 
resolution of the SIS refers to February 1994. 

ellipticals with fewer modeling assumptions. Moreover, the 
calibration status of ASCA has changed significantly over 
the past year and thus the results of previous studies need 
to be verified. We have assembled a sample of 20 early-type 
galaxies (17 E, 3 SO - see Table ^ with high L ^/Lb values as 
listed in Kim, Fabbiano, & Trinchieri ( 1992 ) to insure that 
the X-ray emission is mostly due to hot gas. These galaxies 
have archival ASCA observations of good quality for most of 
the sample. We place a premium on the signal-to-noise ratio, 
S/N, and thus the extraction radii for the spectra are in 
many cases smaller than in previous studies. This attention 
to S/N is vital in order to achieve the global minima with 
XSPEC for our models with many fitting parameters. We fit 
as many temperature components as are demanded by the 
data to achieve a good fit without any restrictions on the 
temperatures. We also examine the reliability of the plasma 
codes in the Fe-L spectral band directly from the sample. 

The organization of the paper is as follows. In §0 we 
describe the details of the observations and data analysis. 
The results of the single-temperatu re fi ts and the empirical 
test of the plasma codes is given in §3.1. We discuss the two- 
temperature (and three-temperature) fits in ^3.2. In §^ we 
discuss the evidence and implications for excess absorption 
in several of the galaxies. Finally, we summarize and give 
our conclusions in 50 . 



2 OBSERVATIONS AND DATA REDUCTION 

The ASCA sateUite (Tanaka, Inoue, & Hoh 1994) consists 
of four detectors each illuminated by a dedicated X-ray tele- 
scope. Two of the detectors are X-ray CCD cameras (Solid- 
State Imaging Spectrometers - SISO and SISl) and two are 
proportional counters (Gas Imaging Spectrometers - GIS2 
and GIS3). Since the X-ray luminosity of the early- type 
galaxies in our sample is thought to be dominated by ther- 
mal emission with T ~ 1 keV n (Forman, Jones, & Tucker 
1985; Canizares, Fabbiano, & Trinchieri 1987; Kim, Fab- 
biano, & Trinchieri 1992) , and the effective area and energy 
resolution of the SIS is substantially greater than that of 
the GIS for energies ~ 0.5 - 1.5 keV (see Table g), the SIS 
data has far better S /N over the important energies for con- 
straining spectral models and for probing the Fe-L emission 
region. 



Throughout the paper we set fcsT — > T. 



That is, the fluxes of bright early-type galaxies above 
2 keV are typically <, 10% of the 0.5-2 keV fluxes (e.g. 
Table ^) which indicates that even though the effective ar- 
eas of the SIS and GIS are comparable for energies ~ 2 — 5 
keV (the SIS being slightly better on average) the important 
constraints on spectral models fltted over ~ 0.5 — 5 keV are 
determined almost completely by the SIS data because of 
its substantial advantage in effective area for energies be- 
low 2 keV and its superior energy resolution over the entire 
ASCA band (e.g. Buote & Canizares 1997). Since the SIS 
data dominates the constraints on spectral models of early- 
type galaxies, any small differences in the model parameters 
obtained from jointly fltting SIS and GIS data of early-type 
galaxies may be due to mismatches in calibration between 
t he two instruments ; e.g. the gain calibration of the GIS 
([desawa et al. 1997). Hence, we focus on the SIS data, al- 
though throughout the paper we illustrate pertinent results 
using the GIS data for representative galaxies. The reduction 
of the SIS and GIS data was performed using the standard 
XSELECT, FTOOLS, and IRAF-PROS software. 

We list details of the observations in Table |^. Because 
the SISO and SISl (and GIS2 and GIS3) have different en- 
ergy responses and gains their data must be fitted sepa- 
rately. The same applies for different observations of the 
same galaxy since the energy response and gain vary sig- 
nificantly with time. As a result, for a galaxy with multi- 
ple observations, if one of the sequences has a substantially 
smaller exposure time than the others, it will not contribute 
significantly to the spectral fits. (And any residual differ- 
ences caused by incorporating the lower S/N observation 
may just be systematics owing to background subtraction 
or calibration.) Moreover, the aperture size used to extract 
the spectrum for the lower S /N sequence will be necessarily 
smaller than the others which could lead to inconsistencies. 
(We describe our determination of aperture sizes below.) 
Hence, for a few galaxies with multiple observations (NGC 
507, NGC 4374, NGC 5044) we ignored these sequences with 
much worse S/N. 

Similarly, for a given observation the SISO data has 
higher S/N than the SISl data because the screening cri- 
teria for events is much more severe for the SISl. If a great 
disparity in S/N existed between the two detectors then we 
analyzed only the high S/N observation (NGC 4649, NGC 
6876). For NGC 507 we excluded the SISl data because of 
anomalies that we discuss in §^ 

For each galaxy we analyzed the data generated by 
the standard processing software that were filtered with the 
default screening criteria (see The ASCA Data Reduction 
Guide 1997). We further screened the data by excluding 
time intervals of high background from analysis of the hght 
curves. For NGC 720 and NGC 1332 we used the events 
s creen ed according to the criteria in Buote & Canizares 
(1997) modified by the light-curve analysis. 

The final filtered events were then extracted from a cir- 
cular region positioned by eye on the galaxy center. A circu- 
lar region was chosen because the standard software which 
incorporates information on effective area and corrects the 
spectra for vignetting (i.e. FTOOLS task ascaarf) currently 
performs best on circular regions. However, for NGC 4406 
we used a square region because the X-ray emission is highly 
asymmetrical which offsets the benefits of a circular region 
in ascaarf. 
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Table 2. Galaxy Properties 



Name 


TvDe 


z 




loETirt /jT3 
olO tj 


^ u 




Group 










orp" ^ ^ 


(km s"^) 


(lO^^cm"^) 




NGC 499 


SO" 


0.01459 


12.64 


43.74 


237 


0.53 


LGG 24 


NGC 507 


LARO 


0.01642 


11.76 


44.19 


366 


0.53 


LGG 24 


NGO 720 


E5 


00572 


11.15 


43.47 


247 


0.14 


LGG 38 


NGC 1332 


S0~ 


0.00508 


11.29 


43.29 


347 


0.22 


LGG 97 


NGC 1399 


EIP 


0.00483 


10.79 


43.46 


310 


0.13 


Fornax 


NGC 1404 


El 


0.00648 


11.06 


43.35 


225 


0.13 


Fornax 


NGC 1407 


EO 


0.00593 


10.93 


43.61 


285 


0.55 


LGG 100 


NGC 3923 


E4-5 


0.00556 


10.79 


43.82 


216 


0.64 


LGG 253 


NGC 4374 


El 


00334 


10.23 


43.67 


287 


0.26 




NGC 4406 


E3 


-0.00076 


10.02 


43.76 


250 


0.26 


V 11 


NGC 4472 


E2 


0.00290 


09.32 


44.04 


287 


0.16 


V 11 g*-" 


NGC 4fi3fi 


E0+ 


u .uuouo 


10.50 


43.58 


191 


0.17 




NGC 4649 


E2 


0.00471 


09.83 


43.83 


341 


0.24 


Virgo 


NGC 5044 


EO 


0.00898 


11.87 


43.75 


234 


0.50 


LGG 338 


NGC 5846 


EO-1 


0.00608 


11.13 


43.77 


278 


0.42 


LGG 393 


NGC 6876 


SBO- 


0.01338 


12.45 


43.66 


231 


0.52 


LGG 432 


NGC 7619 


E2 


0.01269 


12.17 


43.79 


337 


0.50 


LGG 473 


NGC 7626 


EIP 


0.01142 


12.17 


43.70 


234 


0.50 


LGG 473 


IC 1459 


E3 


0.00564 


10.96 


43.53 


308 


0.12 


LGG 466 


IC 4296 


EO 


0.01255 


11.43 


44.12 


323 


0.43 


LGG 353 



Morphological types and redshifts are taken from RC 3. Tot al blue magnitudes (B^^ ) and central 
velocity dispersions (o-q) are taken from Faber et al. (198E) except for NGC 1332 (Dalle Ore et 



al. 1987). Th e lurn inosities (Lb) corresponding to are com puted using the distances from 
' ' Kim (|l 995|) f or Hq = 50 km/s/Mpc. 

LGG refers to Garcia 



Faber et al. (1989) as reported in Eskridge, Fabbiano, 
Galac tic Hydrogen column densities (Afjj) are from Stark et al. (1992) 
(19931). 



Our guiding principle for choosing the radius of the cir- 
cle for a particular galaxy was to maximize the S/N of the 
background-subtracted counts within the ~ 0.5 — 5 keV en- 
ergy band. We modified the S/N criterion according to other 
considerations. We tried to limit as much as possible the cir- 
cular aperture to a single CCD because the responses vary 
from chip to chip. (For the cases where the source counts 
extended over more than one chip we followed the standard 
procedure and averaged the response matrices according to 
the number of source counts on each chip. The gaps between 
chips were excluded.) For the galaxies having a nearby bright 
source (e.g. NGC 1399 and NGC 1404; NGC 499 and NGC 
507) we chose a higher S/N cutoff to reduce contamination 
from the other source. 



from a deep observation of a blank field is used (e.g. Buote 
& Canizares 1997)|T]. 

Finally, the source and background spectra for each 
galaxy were extracted in 512 pulse invariant (PI) bins after 
correcting for temporal SIS gain variations and the degrada- 
tion due to charge-transfer inefficiency. We mention that we 
made the required corrections to fix problems in the response 
matrices generated by the standard reduction software as de- 
scribed on the ASCA Guest Observer Facility WWW pages 
(as of May 1997). Also, the CIS background templates were 
rise-time filtered (i.e. with FTOOL gisclean) to match the 
screening of the galaxy observations. 



3 SPECTRAL FITTING 



The background for each galaxy was chosen from 
source-free regions of the CCDs. Although it is preferred to 
take the background from the same region of the CCDs as 
the source in order to insure the same vignetting correction 
(i.e. using a background determined from a deep exposure 
of a blank field), the variations of the local background due 
to environment and peculiarities for a particular observa- 
tion (e.g. contamination from the bright Earth, the Sun, the 
South Atlantic Anomaly) outweigh the vignetting consider- 
ation for the SIS. However, since the background increases 
mark edly with distance f rom the center of the field of the 
CIS ( Ishisaki et al. 1997 ), we used the blank fields for the 
background estimates of the CIS data. In any event, the 
results using a local background generally differ unimpor- 
tantly from those when a background template constructed 



It has been known since the operation of the Einstein Obser- 
vatory that the X-ray emission of early-type galaxies mostly 
origin ates from a hot interstellar medium (Forman et al. 
197£), with an additional contribution from the integrated 
emission from X-ray binaries which are most important for 
the galaxies with t he smallest ratios of X -ray to blue-band 
optical luminosity (Canizares et al. 1987). The results from 
the Performance- Verification phase of ASCA showed that 
the brightest ellipticals in X-rays have emission lines consis- 
t ent with those ori ginating from a thin hot thermal plasma 
(Awaki et al. 1994). The form of the spectrum arising from 
discrete sources is not well constrained, but is consistent 



^ Typically the fluxes show the largest differences between the 
blank sky templates and a local background. The fitted temper- 
atures and abundances are mostly unaffected. 
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Table 3. ASCA SIS Observation Properties 



Name 


Sequence ^ 


Exposure 


Count Rate 


R 




R 








(^iU 




I^IU 




(arcmin) 


(kpc) 






OliDU 


OlO 1 


OlOU 


OlO 1 




QTQ1 
oioi 


oioU 


oloi 


IMTT^ /inn 


DIUU ( UUU 


ZO.O 


9^^ (<. 
ZO.D 


7 


7 n 

1 .U 


9 O 

z.y 


9 fl 
Z.O 


7/1 7 
( 4. i 


79 9 
( Z.Z 




DvSUzoUUU 


"^A Q 

o4.y 


OO. / 


1 1 

1 1 .0 


Q 1 

y . 1 


O.D 




09 Q 

yz.o 


Q7 a 
O / .O 


IN ou / 


DIUU / UUU 


9"^ 




99 n 
zz.u 




A 9 
4.Z 




191'? 
IZl.O 




■\T^/-i von 


(^nnn/i nnn 
DUUU4UUU 


•^n 9 
ou.z 


"^9 n 


9 ^ 

z.o 


1 A 

1.4: 


9 

Z.O 


9 n 
z.u 


9 "3 7 


ion 
iy .u 


Tvr r^r^ 1 Q'io 


DvSUzZUUU 


UO.U 




1 s 

l.O 


1 n 

l.U 


9 7 
Z. ( 


9 1 
Z. 1 


99 Q 
ZZ.O 


1 7 Q 


IN vjxv^ ioyy 


oUUooUUU 


14 9 


-1-U.O 


34 4 


oo.o 


A O 

4.y 


/I A 
4.4 


oo.y 


04.0 




oUUoyuuu 




10. 1 


OD. Z 


1 Q 

10. y 


/I p. 
4.D 


9 


oD.o 


9f^ Q 
ZO.O 


Tvrr^ 1 /I n/1 


oUUooUUU 


14 9 


1 n 

lU.O 


?i. 9 
o.z 


o. y 


9 7 
Z. 1 


9 

Z.O 


91 zl 
Zl .4 


1 H 9 
lO.Z 




oUUoyuuu 


14:. O 


1 1 

lO. 1 


o.u 


4:.U 


9 

Z.O 


9 n 
z.u 


1 O fl 

iy .O 


10. o 




DOUZiUUU 


•^1 7 


Ol .o 


4 7 

4:. / 


7 
o. / 


O.u 


9 7 
Z. / 


QO Q 

oU.o 


97 Q 


IN VjrL/ oyzo 


ozUiOUUU 


90 ^ 
ZU.O 


19 7 
IZ. I 


1 7 
1. 1 


1 '\ 

l.O 


9 n 
z.u 


9 n 
z.u 


9/1 1 
Z4. 1 


9/1 1 
Z4. 1 


In 4o / 4 


DUUU ( UUU 


-LO.O 


19 7 
IZ. i 


4 n 

4:.U 


4 4 


9 n 
z.u 


9 n 
z.u 


1 7 
iO. t 


1 7 
iO. f 


IN 44U0 


OUUUOUUU i^4ccaj 


1 1 n 

1 l.U 


Q 1 

y . 1 


OO.U 


9fi 9 
ZO. Z 


0.0 


K O 
O.U 


A'^ 9 
40. Z 


oy .O 




DUUUDUUU i^zcca^ 


U.O 


^ 9 
o.z 


ou.z 


90 
zu.o 


O.O 


4.0 


A^ « 
41 .D 


Qfi 1 
oO. i 


i\T r^r^ /I /1 79 


Duuzyuuu 


ID.O 


IZ. 


OD.Z 


9^^ 7 
ZO. / 


A 1 
4. 1 




Q9 9 
OZ.Z 


97 
Z / .O 




OUUoUUUU 


1 7 n 


13 2 


•^9 


90 9. 
ZU.O 


4.3 


2.4 


33.8 


18.8 


IN 4Dv5D 


DUUoZUUU 


33.2 


24.8 


40.6 


27.6 


5.0 


4.6 


39.3 


36.1 


JNCjO 4649 


dUUUdUUU (2ccaj 


16.7 




13.9 




3.0 




23.6 






61005000 (4ccd) 


19.2 




13.8 




3.0 




23.6 




NGC 5044 


80026000 + 80026010 


16.3 


12.2 


117 A 


78.4 


5.5 


5.1 


100.0 


92.3 


NGC 5846 


61012000 


30.3 


16.7 


19.9 


11.5 


4.0 


3.5 


53.3 


46.6 


NGC 6876 


81020000 


18.6 




2.5 




2.0 




53.4 




NGC 7619 


63017000 


52.5 


48.5 


4.7 


3.4 


3.0 


3.0 


54.5 


54.5 


NGC 7626 


63017000 


52.5 


48.5 


0.9 


1.1 


2.0 


2.0 


36.3 


36.3 


IC 1459 


60005000 


23.8 


17.4 


3.9 


1.7 


2.5 


2.0 


23.3 


18.7 


IC 4296 


61006000 


32.9 


31.3 


4.7 


2.7 


3.1 


2.6 


70.7 


59.3 



The exposures include any time filtering. The count rates are background subtracted within the circle 
of radius R. For NGC 4406 we give the half-diagonal of th e rectangular regions. Aperture radii in kpc 
are co mputed using the distances from Faber et al. (1989) as reported in Eskridge, Fabbiano, & Kim 
(|l995|) for Ho = 50 km/s/Mpc. 



with a high-temperature T ^ 5 keV thin thermal plasma 
(Kim et al. 1992; Matsumoto et al. 1997). 

Hence, we concentrate on thin thermal models modified 
by photo-electric absorption arising from the ISM of our 
galaxy. We allow for the p ossibility of intrins ic absorption 
as seen in galaxy clusters (White et al. 1991) if the fit to 
the galaxy spectrum is significantly improved. Since the hot 
gas need not be in a single phase and the X-ray emission 
may contain a significant contribution from discrete sources, 
we fit as many absorbed temperature components as are 
required by the ASCA data. (As stated in §|^, we do not 
assume the presence of a discrete (or other hard) component 
with a predetermined functional form for all galaxies.) 

To model the emission of a thin thermal plasma we 
primarily use the MEKAL code which is a modification of 
the original MEKA code (Mewe, Gronenschild, & van den 
Oord 1985; Kaastra & Mewe 1993) where the Fe-L shell 
transitions crucial to the X-ray emiss ion of ellipticals have 
been re-calculated (Liedahl et. al 1995). We take solar abun- 



dances according to Anders & Grevesse (1989) and photo- 



electric absorption cross sections according to Balucihska- 
Church & McCammon (1992). Although not as accurate 
as the MEKAL model we also e mploy the Raymond- Smith 
code (Raymond & Smith 1977) since it has been used in 
most previous studies. 

We also consider a cooling-flow model based on the 
MEKAL code where the gas cools continuously from some 
upper temperature, Tmax- The emissivity of the cooling gas 



is described in Johnstone et al. (1992). Although this model 
may be more appropriate for the cores of galaxy clusters 
since it does not include potential heat sources which may 
be important in galaxies, it provides a simple and conve- 
nient exploration of the effects of cooling-flow spectra on 
the ASCA galaxy data. The cooling flow spectrum assumes 
that the gas cools at constant pressure. The pressure itself 
does not enter into the calculation, and the same spectrum 
results from different parcels of gas cooling at different pres- 
sures. Provided that the gas at different radii (and thus dif- 
ferent pressures) starts from the same temperature and cools 
before ffowing a significant distance, it is therefore a good 
approximation. In practice, gravitational work will be done 
as the gas ffows, so the true result will be more complicated 
(in a manner which depends on the degree of ffow). An ad- 
ditional temperature component is the simplist modification 
that can be made for this effect (see §3.2.3). 

All of the spectral fi tting was imp lemented with the 
software package XSPEC ( Arnaud 199^] ) using the mini- 
mization method. In order for the weights to be valid for the 
method we regrouped the PI bins such that each group 
contained at least 20 counts. We restricted the fits to ener- 
gies above 0.5 keV (~ 0.55 keV) because of uncertainties in 
the SIS response at lower energies. We also exclude energies 
> 5 keV because the S/N of the bins at these energies is 
generally too low to significantly affect the fits and these 
energies tend to be dominated by uncertainties in the back- 
ground level of the SIS and to calibration errors ( Gendreau 
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& Yaqoob 1997). The normalizations ol the SISO and SISl 
are fitted separately as are the normalization of multiple ob- 
servati ons of the same galaxy. (For the GIS data we analyzed 



energies 0.8 — 9 keV. 



Irwin & Sarazin 1996), NGC 4636 (|Trinchieri et al 1994| ), 
NGC 4649 and NGC 7619 ( [Trinchieri et al. 19961) ,"1TGC 
5044 (David et al. 1994), and several others ( pavis fc Whit" 



For NGC 507 we found that the spectral shapes of the 
SISO and SISl differed substantially. However, we find no 
such discrepancy for the NGC 499 data taken from the NGC 
507 observation. (The N499 and N507 data are take from 
different chips of the SIS.) We compared the NGC 507 data 
with those from a lower S/N observation (61007010) and 
found that the SISO spectra are consistent. (The SISl data 
for 61007010 have too low S/N to make a useful comparison.) 
The GIS data give fitted temperatures and abundances and 
reduced in excellent agreement with the SISO data. Thus 
we restrict our analysis of NGC 507 to the SISO data. (We 
mention that we find no similar SISO/SISl discrepancy for 
the remaining 19 galaxies.) 



3.1 Single-component fits 

3.1.1 Results 

In Table ^ we list the results of simultaneously fitting a 
single absorbed MEKAL model to the available SISO and 
SISl spectra of each galaxy listed in Table ^. In Figure |l| we 
show the spectral fits of NGC 1332 and NGC 4649 (among 
best S/N for lowest L^/Lb) and NGC 4472 and NGC 5044 
(among best S/N for highest Lx/I/b). (For those galaxies 
with multiple observations we show the observation with the 
longest exposure.) In almost all cases the fits are formally 
quite poor with reduced > 1-5. For the X-ray brightest 
galaxies (and those with largest L^/Lb) the residuals of the 
fits tend to be most pronounced around 1 keV, whereas the 
fainter galaxies (and those with smallest L^/Lb) tend to 
show t he largest significant residuals at higher energies. E 
3 keV. I 



1996). We discuss in detail the evidence and implications for 
this excess absorption in ^ 

In terms of both and the fitted parameters the 
Raymond-Smith model gives broadly comparable results to 
the MEKAL model when fitted over the ~ 0.5 — 5 keV range. 
(In agreement with many previous studies we find that the 
Raymond-Smith model tends to give higher temperatures 
and smaller abundances than the MEKAL model.) As ex- 
pected, the models differ in the character of their residuals 
over the Fe-L region ~ 0.7 — 1.2 keV. One effect of these 
differences is that the Raymond-Smith fits are not usually 
improved as much by letting A'h be a free parameter. 

We also display in Table ^ the results for the ab- 
sorbed cooling-flow model. For eleven galaxies the cooling- 
flow model with or without excess absorption fits the spec- 
tra with values comparable to the MEKAL model. For 
the other nine galaxies (IC 1459, IC 4296, NGC 1399, NGC 
1407, NGC 4472, NGC 4649, NGC 5044, NGC 5846, NGC 
7626) the cooling-flow model fits the galaxy spectra better 
than the MEKAL model. In the cases of IC 4296, NGC 4472, 
and NGC 5044 the improvement in the fits is dramatic (see 
Figure ^) . For the galaxies where allowing for excess absorp- 
tion improved the fits to the MEKAL model, the fits with 
the cooling-flow model also improved with excess absorption 
and typically required a larger amount of absorption than 
did the MEKAL model. 

The mass deposition rates of the cooling gas. A/gas , are 
mostly a few solar masses per year consistent with results 
from Einstein (Nulsen, Stewart, & Fabian 1984; Thomas et 
al. 1986). We find generally good agreement with the mass 
deposition rates obtaine d from ROSAT st udies for NGC 507 
(iKim fc Fabbiano 1995t and NGC 1399 (|Rangarajan et al ' 



The temperatures and abundances determined from the 
single-temperature fits are qualitatively consistent with pre- 
vious ASCA studies of the ellipticals in our sample (Awaki 
et al. 1994; Matsushita et al. 1994; Loewenstein et al. 1994; 
Fukazawa et al. 1996; Buote & Canizares 1997; Arimoto et 
al. 1997; Matsumoto et al. 1997; Iwasawa, White, & Fabian 
1997) modulo expected small systematic differences due to 
different source region sizes, background region definitions, 
thermal emission codes used, current status of SIS calibra- 
tion etc. In particular, our single-temperature models re- 
produce the extremely sub-solar abundances found in the 
above studies. All the galaxies have best-fit Z < 0.4:Zq and 
most have Z ^ 0.15Zq. The mean abundance and standard 
deviation are (Z) = 0.19 ± 0.12 Zq. 

For half of the galaxies in our sample the fits improved 
significantly if A^h was allowed to be a free parameter. This 
improvement affects only the energy bins for E ~ 0.5 — 1.0 
keV. In each of these cases, significant excess absorption 
above the Galactic value is indicated. The excess absorption 
we find for these galaxies agrees qualitatively with those 
determined from the previous ASCA studies listed above. 
However, with RO SAT only NGC 1399 sho wed evidence for 
excess absorption ( Rangarajan et al. 1995 ) while several of 
the galaxies which we find to require excess absorption were 
not found to require it wit h ROSAT: NGC 499 and NGC 
507 (|Kim fc Fabbiano 1995|), NGC 4472 (Forman et al. 1993; 



1995), although the Rangarajan et al. value of ~ 1.7Mq 
yr~ for NGC 1399 is about 50 per cent larger than our 
value. We find a large r value of Mgas for NGC 4406 than 
Iwasawa et al. (1997a) because our value refers to a larger 
extraction region. 

For NGC 5044 we obtain a very large Mgas ~ 75Mq 
yr~^ more appropriate for cluste rs of galaxies (e.g . Fabian 
1994). A previous ROSAT study ([Pavid et al. 1994| ) of NGC 
5044 found Afgas « 20 - 25 Mq yr"^ The larger value we 
obtain for NGC 5044 occurs because of the large amount 
of excess absorption required by ASCA whereas for ROSAT 
the absorption is consistent with the Galactic value. 

Interestingly, the abundances obtained from the 
cooling-flow models are generally much larger than found 
from the regular MEKAL models. For the cooling-fiow 
galaxies mentioned above the metallicities all have best-fit 
values Z J; 0.4Zq except for IC 1459. If we consider all of 
the galaxies, the cooling-flow models give a mean and stan- 
dard deviation, (Z) = 0.6 ± 0.5 Zq. 

We mention that the larger values of Tcf with respect 
to those derived from the MEKAL model do not necessar- 
ily imply different gravitating masses. The value of Tcf is 
just the upper temperature from which the gas cools at con- 
stant pressure. The proper temperature to use for determin- 
ing a virial mass is a (gas) mass-weighted temperature. It 
should also be emphasized that there could exist compli- 
cated situations in a multi-phase cooling flow where cold 
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Table 4. Single-Component Fits 
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Results of fitting either one MEKAL model or a cooling-flow model each modified by photo-electric absorption. If 
allowing Nj^ to be free significantly improved the fit, we give the best-fit value of Nji and the change in denoted 
by Ax^(-/Vh). Otherwise "• • •" denote that Nji is kept at its Galactic value. For fits where < 1-5 we also list 90% 
confidence limits on one interesting parameter. The gas mass deposition rates, Mg^s, are quoted in the 0.5-5 keV band 
for the SISO. For galaxies with multiple observations, the highest S/N observation is typically quoted. 
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NGC 4649 NGC 5044 




13 5 12 5 

channel energy (keV) channel energy (keV) 



Figure 1. Absorbed single-temperature (MEKAL) fits to the spectra of NGC 1332, NGC 4472, NGC 5044, and NGC 4649. Only the 
SISO data is shown. 
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gas which has dropped out of the flow is supported by the 
hotter phases (e.g. via magnetic fields). Again, a (gas) mass- 
weighted temperature, or an appropriately adjusted mean 
molecular weight (which will be higher than the single-phase 
hot gas), must be used in the hydrostatic equation. This 
leads to a smaller gravitating mass than that inferred from 
the equation of hydrostatic equilibrium using Tcf. These 
comments als o ap ply to the two-temperature models dis- 
cussed later (§3.2). 



Interpretation of the fitted parameters for IC 1459 must 
be treated with caution. The high temperature, especially 
for the cooling-flow model, implies an implausibly large mass 
for an elliptical of its luminosity. Rather, emission from dis- 
crete sources with a likely contribution from a central AGN 
(due to it being a radio source) suggest that emission from 
hot gas does not dominate the emission. This picture is con- 
sistent with the relatively small value of L^/ Ls of I C 14 59 
and the results of two-temperature fits discussed in §3.2 



Finally, when jointly analyzing the SIS-f GIS data we 
obtained essentially the same answers as when analyzing 
only the SIS data. However, some differences when analyz- 
ing the GIS data alone deserve mention. First, for the high 
Lx/I/B galaxies such as NGC 507 and NGC 4472, the re- 
quired excess absorption depended sensitively on the lower 
energy bound used in the fits. For Eiow = 0.8 keV we ob- 
tained best-fit Nn = 0.01 x 10^^ cm'^ for NGC 507 and 
A^H = 2.3 X 10^^ cm"^ for NGC 4472 whereas for Ei^w = 1 
keV we obtained A"h = 1.8 x 10^^ cm'^ for NGC 507 and 
A^H = 1.6 X 10^^ cm"^ for NGC 4472; i.e. better agreement 
with the SIS data is seen for Eiow = 1 keV. Second, for the 
low Lx/I/B galaxies such as NGC 1332 and NGC 3923 we 
find that the best-fitted temperatures are approximately 2 
keV as opposed to about 0.6 keV with the SIS. These val- 
ues are formally inconsistent within their errors and simply 
reflect the different energy resolutions and different effective 
areas as a function of energy of the SIS and GIS as discussed 
in U 



3.1.2 Assessing the validity of measured abundances 

The very sub-solar abundances obtained from the single- 
temperature MEKAL and Raymond-Smith fits has cast 
doubt on the reliability of the plasma codes, in particular 
the accuracy to which the emission in the Fe-L spectral re- 
gion, 0.7 — 1.4 keV, is modeled (e.g. Arimoto et al. 1997). 
We have empirically tested the reliability of the MEKAL 
code by examining the behavior of the fitted abundances 
and temperatures when the Fe-L spectral region is excluded 
from the fits to the SIS data. This procedure is motivated 
by the fact that for gas with T ~ 1 keV the peak in the 
galaxy spectrum occurs in the Fe-L region: the location of 
this peak is one measure of the temperature. Outside the 
Fe-L region the continuum level provides an independent 
measure of the temperature. Comparing the temperatures 
(and abundances) determined inside and outside the Fe-L 
region directly tests the consistency of the plasma code. 

We analyzed the restricted energy range, E = 0.5 — 2 
keV, to limit the effects of a possible hotter component in the 
hot gas and/or a hard component due to discrete sources. 
We fitted the MEKAL model with free A'h to the following 
two cases: (in) - the Fe-L region and (out) - the energies 
excluding the Fe-L region over 0.5-2 keV. For both (in) and 
(out) the initial parameters in the fits were specified by the 
model obtained from fitting the entire 0.5-2 keV range. For 
these fits we defined the Fe-L region to be i? = 0.7—1.2 keV. 
Although Fe-L contributes up to energies ~ 1.4 keV, the 
emission there is comparable to those from other elements 
(most notably Mg). 

Let us suppose that the MEKAL model produces too 
much Fe-L emission for a given T and Z^. The fits can sup- 
press the excess Fe-L emission by either increasing T or de- 
creasing Z or both. Hence, in this case we would expect fits 
to region (in) to have higher T and/or smaller Z than fits to 
(out). We would expect the opposite behavior if the MEKAL 

■t- Relative abundances are fixed at their solar values. 
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model produces too little Fe-L emission for a given T and Z. 
(These arguments assume a single temperature component 
dominates the emission over 0.5-2 keV.) 

In Figure ^ we plot temperatures and abundances for re- 
gion (out) against those of region (in). For all but two galax- 
ies (IC 1459, IC 4296) the temperatures show little scatter 
along the line Tin = Tout giving a strong indication of the 
validity of the MEKAL code in the Fe-L region. The high 
inferred temperature of IC 1459 (see Table 0) accounts for 
its large value of Tout while the significant contribution from 
a higher temperature component ov er re gion (out) likely ex- 
plains the position of IC 4296 (see ^3.2). 

The abundances show considerable scatter, but there 
appears to be an excess of galaxies above the line Zin = 
^out; i.e. Zout tends to exceed Zin in our sample, though the 
effect is not large. Given the consistency of the temperatures, 
this excess Zout may indicate that the abundances of the 
elements dominating the line emission of region (out) (i.e. 
Si) are slightly over-abundant with respect to Fe. Once again 
IC 1459 is the outlier (leftmost point). 

Finally, in Figure ^ we plot the ratio of temperature 
to abundance for region (out) versus region (in). If the Fe-L 
emission of the MEKAL model is incorrect, then T/Z should 
tend to be larger in region (in) if MEKAL predicts too much 
Fe-L emission or smaller if MEKAL predicts to little Fe-L 
emission; i.e. the slope in the T/Z (in) - T/Z (out) plane will 
be negative if the MEKAL model incorrectly predicts the 
Fe-L emission. However, no such trend is found in Figure ^ 
Rather, T/Z for the regions tend to track each other, albeit 
with large scatter which is due mostly to the scatter in the 
abundances. 

Hence, over 0.5-2 keV the MEKAL code gives very con- 
sistent results when applied to only the Fe-L region (in) and 
to the region excluding the Fe-L lines (out). The very sub- 
solar abundances determined from the single-temperature 
fits of the MEKAL model are thus not the result of er- 
rors in the Fe-L portion of the plasma code. As we show 
in the next section, these sub-solar abundances arise from 
fitting a single-temperature model to an intrinsically multi- 
temperature spectrum. 

(We have performed the same tests using the Raymond- 
Smith model. The results are qualitatively similar to those 
found with the MEKAL model, although the scatter in 
the parameter relationships in regions (in) and (out) with 
the Raymond-Smith model is larger. Also Tin tends to ex- 
ceed Tout for the Raymond-Smith model. Although the 
Raymond-Smith model, as expected, gives more evidence 
for problems in the Fe-L regions than the MEKAL model, 
the differences do not appear to be huge. This is sup- 
ported by the qualitative agreement that we obtained for 
the temperatures and abundances of the single-temperature 
MEKAL and Raymond-Smith models, though the higher 
temperatures and lower abundances of the Raymond-Smith 
model found in the previous section are consistent with over- 
prediction of the Fe-L emission. As we discuss below, a 
large discrepancy between the MEKAL and Raymond-Smith 
models is not found for the two-temperature fits as well.) 



3.2 Two-component fits 

Since most of the single-temperature fits with the MEKAL 
model clearly provide unacceptable fits to the ASCA spec- 




T. /Z. 

in' in 

Figure 4. As Figure ^except that the ratios of temperature to 
abundance are plotted. 



tra, we examined whether the fits would improve signifi- 
cantly upon adding another temperature component. Our 
procedure is to begin with the best-fitted single-temperature 
model with Galactic A^h and then add another absorbed 
temperature component with A^h and Z tied to the first 
component. The normalization of the second component is 
initially set to zero. After fitting this two-temperature model 
we then allowed A^h to be free to see if the fit could be fur- 
ther improved. 

The process of fitting these two component models to 
the ASCA spectra required some care to insure that the 
global minimum was achieved by XSPEC for the galax- 
ies having lowest S/N in our sample; i.e. for approximately 
half of the galaxies. Typically, upon initially fitting the two- 
temperature model to these lower S/N galaxies one finds 
that the first component, let us call it the cold compo- 
nent with temperature Tc, has best- fit Tc 1 keV simi- 
lar to what was obtained with the single-temperature fits. 
The temperature of the second component, let us call it the 
hot component with temperature Th, typically has best-fit 
Th 10 keV. Also, the abundances tended to differ insignif- 
icantly from their single-temperature values. 

We found, however, that for many of these galaxies 
these fits were not the global minima. By taking these best- 
fit results and then initially setting Th to some relatively 
small value (e.g. 1 — 2 keV) and the abundances to some 
relatively large value Z ~ (1 — 2)Zq and then fitting again 
we found that a deeper minimum was usually obtained. The 
results of this new deeper minimum were generally smaller 
values of Th and larger values of Z, though in most cases 
the 90 per cent confidence regions overlapped with those of 
the previous minimum. 

For example, the initial fit for IC 4296 yielded (Tc = 
0.81 keV, Th = 18.4 keV, Z = 0.20Zq, = 93.8). By 
then resetting Th and Z to small and large values respec- 
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Table 5. Two- Temperature Fits 



Name 


Nn 




Tc 


Th 


z 


/c//h 


Ax'(JVh) 




dof 


y2 




(1021 cm-2) 




(keV) 


(keV) 


(%) 


(0.5-2 keV) 


(0.5-5 keV) 










NGC 499 


2 

^■^-0.4 





fic+O.OS 
°°-0.02 


5.7(> 3.1) 


0.67tO-j^ 


11^ 

'-'--4 


4 


32.2 


231.9 


200 


1.16 


NGC 507 


1 O+0.5 

^^■■'-o.s 





oq+0.25 
"^-0.22 


-^■^"-0.17 


"■"•^-0.20 


n Q9+0.31 

"■'-'^-0.23 


9C+0.32 


10.4 


60.8 


72 


0.85 


NGC 720 







c;7+0.04 
""-0.05 


7.4(> 3.6) 


0.4(> 0.1) 


9 q+1'7 


1 4+0.7 




42.6 


45 


0.95 


NGC 1332 







C.+0.04 
^^-0.05 


14(> 7) 


0.8(> 0.2) 


9 99+0.78 
^•^^-0.22 


"■^"-0.20 




66.9 


63 


1.06 


NGC 1399 


9+"-^ 





yo+0.05 


^■"^-0.09 


1 19+0.43 
-^■-^^-0.23 


"■^"-0.03 


SI +0-03 
"■■^^-0.03 


28.8 


310.9 


251 


1.24 


NGC 1404 







gQ+0.02 
°^-0.02 


9 1 q+0.92 
^- J-'J-0.40 


7+2 
"■ '-0.3 


Q q+1-2 
■^-^-1.3 






168.8 


123 


1.37 


NGC 1407 







72+0.04 
'^-0.04 


3.9(> 2.9) 


n c:o+0.98 
"■^•^-0.36 


^■"^-0.33 


0.87^0;?° 




84.5 


91 


0.93 


NGC 3923 







cc+0.08 
'^^-0.13 


4.2(> 2.2) 


2.0(> 0.1) 


1 sq+0.93 


QQ+0-65 
U.»»_Q 23 




21.8 


20 


1.09 


NGC 4374 







gg+0.05 
""-0.04 


5.0(> 2.1) 


"■^'-0.09 


7 q+19 
'■'^-3.9 


q+10 
•^■^-1.8 




59.3 


49 


1.21 


NGC 4406 


1 0+°-^ 

^■^-0.4 





7n+0.03 
'"-0.07 


1 99+0.30 
-^■^^-0.28 


n c:9+0.22 
"■^^-0.18 


2 q+i-'^ 


2 4+1-2 
^■^-1.0 


11.5 


179.3 


183 


0.98 


NGC 4472 


1 9+0.3 
-^■^-0.3 





'"-0.03 


1 48+0-11 
-^■^''-0.09 


1 10+0.39 


1 +0.11 
^■'-'^-0.08 




56.6 


260.0 


235 


1.11 


NGC 4636 


1 9+0.3 
-^■^-0.2 





oc:+0.02 
"°-0.01 


> 5.3 








46.8 


206.0 


135 


1.53 


NGC 4649 


1 t;+0.6 





7c:+0.04 
'^-0.04 


1 90+°-*2 






1 36+0-17 


21.0 


131.8 


96 


1.37 


NGC 5044 


1 5+'''2 





7n+0.02 
'"-0.03 


1 20+°-*^^ 
-^■^"-0.05 


62+" " 
"■"^-0.08 


n sc;+0-09 
^-^^^-0.08 


r, 79+0.08 
"■ '^-0.06 


82.0 


193.9 


158 


1.23 


NGC 5846 


1 q+0.5 

^■■'-0.4 





cq+0.03 
"•^-0.04 


1 O9+0.18 
-^■'-'^-0.21 


"■^^-0.19 


9 + 0.9 


9 9+0.6 


13.6 


92.8 


104 


0.89 


NGC 6876 







90+014 
^"-0.14 


3.0(> 1.5) 


97+0.66 
"■^'-0.17 


3 5+29 






21.0 


16 


1.31 


NGC 7619 


o 4+1.0 





fic+O.OS 
"""-0.05 


3.1(> 1.8) 


1.2(> 0.3) 


q 3+16.1 


43+3.5 


31.3 


82.6 


95 


0.87 


NGC 7626 







gq+0.07 
"^-0.06 


o q + 1.6 

■'•■^-i.o 


2.0(> 0.4) 


^-"-0.4 


1 l+°-2 




52.8 


49 


1.08 


IC 1459 







gg+0.09 
""-0.10 




"■l^*-0.09 


"-■^"-0.01 


99+0.26 
"■^^-0.03 




58.9 


52 


1.13 


IC 4296 







7cr+0.05 
'■^-0.04 


5.0+1:3 


2.8(> 0.9) 


94+0-47 

"■^^-0.13 


49+0-32 

"■^^-0.05 




86.7 


88 


0.99 



Results of fitting 2 MEKAL models each modified by the same photo-electric absorption. If allowing Nn to be free significantly improved 
the fit, we give the best-fit value of TVh and the change in denoted by Ax2(]Vh). The "■ ■ •" indicate that letting A^h be free did not 
significantly improve the fits. The abundances, Z, are tied together in the fits. The best-fit values and their 90% confidence limits on 
one interesting parameter are given. /c//h is the flux of the absorbed cold component divided by that of the absorbed hot component. 



Table 6. Three- Temperature Fits 



Name 


Ah 


Tl 


T2 


T3 






(10-21 cm-2) 


(keV) 


(keV) 


(keV) 


(Z 


NGC 4472 


1 q+0.3 

1-^-0.3 


q 70 + 0.03 

"•'■^-0.03 


1 34+0.10 
-^"^^-0.08 


> 6 


1.19 


NGC 4636 


1 4+0-5 


53+0-10 


r, 70+0.20 

"■ ''-'-0.07 


8(> 3) 


0.47 


NGC 4649 


^•"-0.7 


"•"'-0.08 


1 15+0-21 


> 5 


0.88 



-0.46 
-0.28 
+0.17 
0.12 
+0.41 
0.33 



/l :/2 


:/3 




dof 


-^red 


(0.5-2 keV) 


(0.5-5 keV) 








21 : 17 : 1 


7.0 : 7.1 ; 1 


246.4 


230 


1.07 


15 ; 12 : 1 


6.1 : 5.0 ; 1 


181.7 


132 


1.38 


9.7 ; 7.3 : 1 


3.1 : 2.7 ; 1 


110.0 


93 


1.18 



Results of fitting 3 MEKAL models each modified by the same photo-electric absorption for three galaxies where the fit is at 
least marginally improved. The abundances, Z, are tied together in the fits. The best-flt values and their 90% confldence limits on 
one interesting parameter are given for the column densities, temperatures and abundances. The best-fltting ratio of the fluxes of 
component 1 to component 3 and the ratio of component 2 to component 3 are given by /i : f2 '■ fs- 



tively, we obtained (after performing this procedure twice) 
best fit values (Tc = 0.75 keV, Th = 5.0 keV, Z = 2.8%, 
~ 86.7) as given in Table h the change in is greater 
than 2.71. The other galaxies showed a smaller effect. Ini- 
tial fits to NGC 720 gave (Tc = 0.57 keV, Th = 41.4 keV, 
Z = 0.18%, = 43.7). The resuhs listed in Table | show 
that the change in x^ is only 1.1, but the best-fit parame- 
ters and their confidence regions (though overlapping) are 
qualitatively different; i.e. the abundance in the deeper min- 



imum is unbounded from above unlike that in the shallower 
minimum. 

In Table ^ we list the results of fitting two absorbed 
MEKAL models to each galaxy in the sample. In Figure ^ 
we plot the two-temperature fits corresponding to the same 
galaxies in Figure |l| In all but a few cases the improvement 
in the fits is substantial and the values of reduced x^ are 
brought to approximately 1.0. (We found that untying the 
abundances and /or absorption of the components did not 
significantly improve the fits.) 
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As expected, the galaxies like NGC 507, NGC 4374, 
and NGC 4406, which already had reduced « 1.0 from 
fits to a single MEKAL model, generally showed the small- 
est improvement when another component was added. For 
NGC 1404 the poor fit is mostly due to residuals near 1.8 
keV in the observation sequence 80039000. The origin of 
these residuals is unclear (they are not seen in the spectrum 
of NGC 1399 from the same sequence), but if we exclude 
this sequence we obtain a reduced ^ 1-0 for the single- 
temperature model and thus, like NGC 507, NGC 4374, and 
NGC 4406, a single temperature appears to produce the 
emission without great (formal) need for another compo- 
nent. 

For NGC 6876, which has the poorest S/N of the 
sample, letting A'h be free substantially improved the fit, 
~ 16.5, but the resulting fitted parameters were very 
poorly constrained. (Because of the low S/N we were unable 
to obtain a definite best-fit result because the x^ space con- 
tained many small isolated minima of similar depth.) The 90 
per cent confidence regions, however, include the parameter 
space obtained when A'h is held at the galactic value, and 
thus we have presented those results in Table pi 



As expected (see when fitting the SIS-I-GIS data 
we found, as with the single-temperature fits of the previ- 
ous section, that the results differed negligibly from the SIS 
results alone. We shall henceforth focus on the SIS results. 

3. 2. 1 Temperatures 

The effect of adding the second temperature component can 
be divided roughly into two classes. For the galaxies with the 
best S/N and typically the largest Lx/ie the two tempera- 
ture components essentially bracket the single-temperature 
value with Tc ~ 0.5 - 1 keV and Th ~ 1 - 2 keV. These 
temperatures are consistent with emission from hot gas. For 
the galaxies having large S/N and large Lx/^b but Th > 2 
keV (NGC 499, NGC 4636, NGC 7619), the cold compo- 
nent dominates the X-ray emission. Hence, the high S/N, 
high Ly^/Ls galaxies either have no evidence for emission 
from a hard component with Th ^ 5 keV or allow for only 
a small fraction (;$ 10%) of the total flux to reside in such 
a component. Rather, multiple phases in the hot gas are 
indicated by the two-temperature (and three-temperature, 
see below) fits. The fits to most of these galaxies are also 
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significantly improved by allowing for excess absorption. In 
g4l we discuss in detail the evidence and implications for 
this excess absorbing material. (We reiterate that the fits 
are not significantly improved by allowing for individually 
fitted abundances for the cold and hot components, or by 
letting the a-process elements be fitted separately from Fe.) 

The other galaxies, those with lowest S/N and small- 
est Lx/I/B, typically have Tc ~ 0.5 — 1 keV (similar to the 
single-temperature value) and Th ~ 5 — 10 keV. The cold 
component is consistent with emission from hot gas while the 
hotter component has temperatures consistent with that ex- 
pected from discrete sources (e.g. Kim et al. 1992); i.e. this 
putative discrete component is found to be most important 
in galaxies that tend to have the lowest L^/ Lb ratios in our 
sample as expected (Canizares et al. 1987). Typically the 



emission from the cold component exceeds that of the hot 
component in the soft energy range, 0.5-2 keV, while the 
emission from each component is roughly equal over the full 
0.5-5 keV band. IC 1459 is the exception as its hot com- 
ponent dominates which again suggests that its emission is 
qualitatively different from others in our sample and may 
indicate the emission from the AGN. 

For the galaxies with Th ~ 3—5 keV having large L^/Lb 
but average S/N in our sample, it is possible that higher S/N 
data will split Th into another soft component and a hard 
component. (This may even happen for the lowest L^/Lb 
galaxies.) Interpretation of the nature of the hotter compo- 
nent in these galaxies must await better quality data. 

NGC 4636 has the worst formal fit of the galaxies in our 
sample which does not have any obvious anomalous features 
in its spectrum (see Figure ^) . Untying the absorptions and 
abundances of the components did not improve the fits. We 
tried grouping the a-process elements (O, Ne, Mg, Si, and 
S) separately from Fe and the rest to no avail. The best im- 
provement in was obtained by adding a third temperature 
component (see Table ^ and Figure ^). 

The fit was not improved greatly and so there may be 
something special to NGC 4636 inhibiting a good fit. A new 
200ks observation of NGC 4636 shows the same anomalous 
spectrum indicating that the features are real and not the re- 
sult of a spurious observation (Buote, Fabian, & Canizares 
1997, in preparation). The addition of the third tempera- 
ture component split Tc into two roughly equal components 
with temperatures between 0.5-1 keV. The flux of the high- 
temperature component remains dwarfed by the soft com- 
ponents. The substantial residuals for E < 1 keV and the 
fact that NGC 4636 has a l ow va lue of Tc for its X-ray lu- 
minosity (see Figure ^ in §3.2.4) suggests the presence of 
additional cold phases. 

Essentially all of the other galaxies are not affected sig- 
nificantly by adding another temperature component, al- 
though NGC 4472 and especially NGC 4649 show some im- 
provement (see Table |^ . In these two cases the added com- 
ponent has a temperature with no upper bound and may 
indicate the emission from discrete sources in those galaxies. 
Again, the softer components dominate the hard component. 

Fitting two Raymond-Smith models to the galaxy spec- 
tra gives values that are qualitatively similar to those ob- 
tained from the MEKAL models in many cases, but there are 
notable differences. Generally, the reduced values of the 
Raymond-Smith fits are larger, (x?ed)^^ ^ iXr 



NGC 4472 (x?ed = 1-47) are noticeably poorer than the 
MEKAL fits while NGC 4636 (xL = 2.65) and NGC 4472 
(Xrod = 1.96) have dramatically worse fits. These galaxies 
coincidentally have the best S/N (and among the largest 
Lx/Lb) in the sample. 

When A^H was allowed to be free the fits with the 
Raymond-Smith model were not improved as much as were 
fits with the MEKAL model. Only NGC 507, NGC 4649, 
NGC 5044, NGC 6876, and NGC 7619 showed significant 
improvement but with a lesser change in x^ than found with 
the MEKAL model. 

The fits to the galaxies using the Raymond-Smith code 
also gave qualitatively comparable temperatures to those de- 
termined with MEKAL, although both Tc and Th tend to be 
larger for Raymond- Smfth|5[ Typically, T^^ ~ T^^^^^+0.2. 
For galaxies where Th ~ 5 keV as determined from the 
MEKAL model, the Raymond-Smith fits tend to give Th > 
10 keV; for many of these only lower limits could be found. 
We found that the Raymond-Smith fits gave Th < 2 keV for 
only NGC 1404, NGC 4472, NGC 5044, and NGC 5846. In- 
terestingly, when A^H is allowed to be free for NGC 4472 the 
best-fit value of Th increases without bound, although the 
quality of the fit is not significantly improved. (The abun- 
dances in these two cases are quite different as we mention 
below.) Hence, the nature of the emission of the second tem- 
perature component, be it due to another phase of the hot 
gas or to discrete sources, is dependent on the plasma code; 
i.e. MEKAL tends to give smaller Th and thus more indi- 
cation of additional phases of the hot gas than Raymond- 
Smith. 

Ten of the galaxies in our sample have been previously 
analyzed by Matsumoto et al. (1997) using ASCA data (see 
. These authors fit a two-component model to the SIS and 
CIS spectra where each component is modified by the same 
variable absorption, and the abundances of each component 
are tied together. For IC 1459, IC 4296, NGC 499, NGC 720, 
NGC 4374, and NGC 4636, where we found Tc ~ 0.5- 1 keV 
and Th 5 keV, Matsumoto et al. obtained Tc and reduced 
X^ values (and A^h) in good agreement with our results - 
whether we used two MEKAL models or Raymond- Smith 
models. (We have different reduced x^ values for NGC 507 
because we excluded its anomalous SISl data - see begin- 
ning of ^.) For NGC 4406, NGC 4472, and NGC 4649 Mat- 
sumoto et al. obtained comparable quality fits to us without 
requiring a second hot phase. As mentioned above, this is 
due to fitting Raymond-Smith models instead of MEKAL 
models. 

We also obtaine d ess entially the same Tc and Th of 
Buote & Canizares (1997) who fitted two Raymond-Smith 



components to NGC 720 and NGC 1332. The unbounded 
Th (and low abundance) obtained by Buote & Canizares 
using two Raymond- Smith models agrees with our findings 
using the same models. Our Raymond-Smith fits for NGC 
4374, NGC 4406, NGC 4472, and NGC 4636 give Tc and 
Th in good agreement with the results of Matsushita et al. 
(1994) who fitted two Raymond-Smith components to the 



0.2. However, the fits to NGC 1399 



(Xrcd 



1.43) 



and 



CIS data. 

Our results for NGC 5044 appear to be consistent 



3 This is consistent with the Raymond-Smith model p redict ing 
too much Fe-L emission for a particular T and Z - see §|3.1.4 
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with those from Fukazawa et aL (L996) who fitted single- 



temperature models to ASCA SIS and GIS data in the re- 
gions R < 2' and A' < R < 6' and found x?cd > 1-42 for 
a suite of 4 plasma codes. Since the residuals of the fit of 
their single-temperature Raymond-Smith model (see their 
Figure 4) are very consistent with our findings (e.g. with 
the MEKAL model in Figure a two-temperature model 
with both temperatures less than 2 keV would appear to be 
required for an acceptable fit to their spectra. 



3.2.2 Abundances 

The metal abundances obtained from the two-temperature 
(and three-temperature) MEKAL models are substantially 
larger than those obtained from the single-temperature mod- 
els. For 18 out of 20 galaxies we find that the best-fit 
metal abundances for the two-temperature models have 
Z > 0.4Zq. The mean and standard deviation of the abun- 
dances for the entire sample are {Z) — 0.9 ± 0.7 Zq. The 
two galaxies (IC 1459 and NGC 4374) having very subsolar 



abundances also have among the lowest S/N spectra in our 
sample]^. 

This systematic increase of Z derived for the two- 
temperature models with respect to single-temperature 
models occurs because of an intrinsic bias specific to X- 
r ay sp ectral fitting. As noted in §2.2 of Buote & Canizares 
(Fggi) regarding the analysis of the ROSAT PSPC spec- 



trum of NGC 720 with the Raymond-Smith code, when a 
galaxy spectrum which is intrinsically described by a two- 
temperature solar abundance model is fitted with a single- 
temperature model, the fitted abundance will necessarily 
have a value that is very subsolar; i.e. the very subsolar 
abundances are an artificial systematic bias of fitting single- 
temperature models to intrinsically multi-temperature spec- 



^ We note that it would be better to quote a mean abundance 
weighted by the uncertainties in the individual measurements. 
However, the error bars are highly asymmetric and the sizes of 
the errors are highly correlated with the measured value of Z. 
Thus, a simple weighting by standard errors is inappropriate. 
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tra. We verified this result using simulated two-temperature 
ASCA spectra with MEKAL models for the highest S/N 

galaxies in our present sample|U]. (We mention that for the 
galaxies where excess absorption improved the fits, when the 
absorption is fixed to its Galactic value the abundances are 
usually about the same as for the fits allowing for excess 
absorption, though in some cases the abundances are con- 
siderably larger; e.g. for NGC 4649 we obtain Z = 0.89Zq 
when allowing for excess absorption but Z = 1.24^0 for 
Galactic A^h. Thus the relatively large abundances are not 
the result of the (in some cases) large required Nh suppress- 
ing the excess line emission at lower energies.) 

Hence, the very subsolar abundances obtained by pre- 
vious studies with ROSAT (e.g. Davis & White 1996) 
and ASCA (e.g. Arimoto et al. 1997) which fitted single- 
temperature models to the X-ray spectra of ellipticals ap- 
pear to be a fitting artifact. However, even recent studies 
with ASCA which fitted two-temperature models obtained 
typical abundances that are signifi cantly less than solar. For 
example, Matsushita et al. (1994) analyzed a small sample 
of early-type galaxies and obtained Z ~ OAZq for NGC 
4406, NGC 4472, NGC 4636, and Z ^ 0.1% for NGC 
4374 by fitting two Raymond-Smith components to CIS 
data. Our Raymond-Smith results produce similar values 
except for NGC 4472 where we obtain Z = 1.3Zq. We ob- 
tain qualitatively similar results for these galaxies with fits 
to two MEKAL models. If we let Ah be free for the fit to 
NGC 4472 with two Raymond-Smith models, then we ob- 
tain Z = 0.5Zq in much better agreement with Matsushita 
et al. (although the fit is not improved significantly). 

These four galaxies, however, have somewhat smaller 
abundances than our entire sample. For the fits with two 
Raymond-Smith models we obtain for the whole sample an 
abundance with mean and standard deviation, (Z) = 0.7 ± 
0.6 Zq. This value is systematically smaller than obtained 
from the fits with two MEKAL models, but the relatively 
s mall systematic difference is consistent with our findings in 



§3.1.2 



However, these abundances are still systematically 
larger than the abun dances {Z ^ 0.4Zq) obtained by Mat- 
sumoto et al. (19971) who analyzed a sample of 12 early- 
type galaxies (see fjll). The disagreement with our results 
probably is due to the difi'erent modeling procedures. That 
is, although Matsumoto et al. use a variable temperature 
Raymond-Smith model for their soft component, they use 
a fixed, 12 keV, Bremsstrahlung temperature for the hard 
component. (And they let Ah be free in all cases, regard- 
less of whether the fit demands that it be different from 
the Galactic values.) We explained near the beginning of 



§3.2 that larger values of Th tend to imply smaller values of 



Z, especially for the galaxies having spectra with the lowest 
S/N. Considering that Matsumoto et al. generally use larger 
extraction apertures than us (particularly for the SISl) their 
spectra have lower S/N and thus these effects become more 
pronounced. 



II Some ROSAT studies cautioned that poor quality spectra sim- 
ply could not distinguish between single-temperature models with 
very sub-solar abundances and two-temperature models with so- 
lar abundances (e.g. Trinchieri et al. 1994) owing to the larger 
number of free parameters in the latter case. 



Thus, the modeling differences and lower S/N data 
likely account for the differences between our present study 
and those of Matsumoto et al. (The updated ASCA calibra- 
tion used in our paper also contributes to these differences.) 
Since the MEKAL model is more accurate than Raymond- 
Smith, and we believe there is no reason to bias the fits with 
a preconceived hard component required to be present in 
the same form in all galaxies, and our data have generally 
higher S/N, the approximately solar abundances found in 
this paper would appear to be favored. Possible abundance 
gradients do not resolve our discrepancy with previous stud- 
ies because we obtain approximately solar abundances even 
for some of the large S/N galaxies in our sample which have 
large extraction radu (e.g. NGC 1399, NGC 4472). Approx- 
imately solar abundances are consistent with multi-phase 
models for the evolution of the hot gas in ellipticals (e.g. 
Fujita, Fukumoto, & Okoshi 1997) but are still inconsis- 
tent with the highly super-solar abundances predicted by 
standard enrichment models (e.g, Ciotti et al. 1991; also see 
Arimoto et al. 1997). 



3.2.3 Cooling flows 

We also fitted the galaxy spectra to a two-component model 
consisting of a cooling-flow model for one component and 
a MEKAL model for the other. The absorption and abun- 
dances were handled as with the two MEKAL models above. 
In all cases we found that these models gave fits of essentially 
the same quality as obtained with the 2 MEKAL models. 
Generally, the fitted temperatures and abundances behaved 
similarly as well. For galaxies with both temperature com- 
ponents less than about 2 keV the MEKAL component can 
be thought to represent either the ambient gas phase or the 
gravitational work done on the flow, and the cooling-flow 
component the distributed mass drop-out term. It should 
be emphasized that Tcf is the maximum temperature from 
which the gas cools (i.e. not the emission-weighted temper- 
ature). 

For example, we found for NGC 5044 the temperature of 
the cooling-flow component to be, Tcf = 1.6 keV, that of the 
MEKAL component to be, Tm = 0.94 keV, and the abun- 
dance to he Z = 0.94Zq; as with the single-temperature 
flts, the abundances obtained with cooling-fiow models are 
larger. The fiux of the cooling-flow component is dominant, 
/cf//a/ ~ 3.6 (0.5-2 keV). The inferred mass deposition 
rate, Mgas ~ 47Mq yr~^, is less than that obtained in the 
single-component case since the flux is now shared with the 
other component. 



3.2.4 Luminosities 

In Table ^ we give the unabsorbed luminosities of the cold 
and hot components obtained from the flts of the two 
MEKAL models to the galaxy spectra as well as the ra- 
tio of total X-ray luminosity, Lx = + Lx , to Lb. (We 
give only the best-fit values of Lx since the statistical un- 
certainties are less than a few percent. The statistical errors 
on the relative contribution of L^ and L^ are very nearly 
those of /c//h in Table |.) Because of ourprocedure of max- 
imizing the S/N of the extraction radii (§0) our luminosities 
tend to be somewhat smaller than found in previous studies 
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Table 7. Luminosities 
















Name 


log 


10 




logic L^/L 


B 






(0.5-2 keV) 


(0.5-5 keV) (0.5-2 keV) 


(0.5-5 keV) 


(0.5-2 keV) (0.5-5 keV) 




NGC 499 


42.55 


42.57 


41.39 




41.70 


-1.16 


-1.12 




NGC 507 


42.19 


42.21 


42.66 




42.76 


-1.40 


-1.32 




NGC 720 


40.76 


40.77 


40.30 




40.58 


-2.58 


-2.48 




NGC 1332 


40.43 


40.44 


40.07 




40.43 


-2.70 


-2.55 




NGC 1399 


41.30 


41.32 


41.66 




41.77 


-1.65 


-1.56 




NGC 1404 


41.23 


41.24 


40.64 




40.80 


-2.02 


-1.97 




NGC 1407 


40.97 


40.99 


40.75 




41.01 


-2.43 


-2.31 




NGC 3923 


40.89 


40.90 


40.58 




40.83 


-2.76 


-2.65 




NGC 4374 


41.04 


41.06 


40.17 




40.45 


-2.57 


-2.51 




NGC 4406 


41.78 


41.79 


41.29 




41.37 


-1.86 


-1.83 




NGC 4472 


41.59 


41.61 


41.43 




41.54 


-2.22 


-2.17 




NGC 4636 


41.97 


41.98 


40.28 




40.69 


-1.60 


-1.58 




NGC 4649 


41.41 


41.42 


41.06 




41.20 


-2.26 


-2.20 




NGC 5044 


42.80 


42.81 


42.82 




42.89 


-0.64 


-0.60 




NGC 5846 


41.96 


41.97 


41.46 




41.54 


-1.69 


-1.66 




NGC 6876 


41.51 


41.55 


41.00 




41.23 


-2.03 


-1.94 




NGC 7619 


42.25 


42.26 


41.09 




41.31 


-1.51 


-1.49 




NGC 7626 


41.06 


41.07 


40.82 




41.05 


-2.44 


-2.34 




IC 1459 


40.47 


40.49 


40.71 




41.03 


-2.62 


-2.39 




IC 4296 


41.53 


42.00 


41.54 




41.81 


-2.29 


-1.90 




Intrinsic (unabsorbed) luminosities of the cold component (i£^) and the hot component {L^) in 




erg s corresponding to the two-temperature models of 


I'able bI We used the distances from Faber 




et al. (|l98q) 


as reported in 


Eskridge, Fabbiano, & Kim ( 


1995 


k for Hq = 


50 km/s/Mpc. 
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Figure 7. (left) The total 0.5-5 keV band X-ray luminosity versus Lb and (right) 0.5-5 keV band luminosities of the cold components 
(boxes) and hard components (crosses). 



which is most significant for galaxies having the lowest S /N 
observations. However, the 0.5-5 keV fluxes determined in 
this paper differ by ^ 25% for the corresponding galaxies 
studied by other authors using larger aperture sizes (Awaki 
et al. 1994; Matsushita et al. 1994; Loewenstein et al. 1994; 
Buote & Canizares 1997; Arimoto et al. 1997; Matsumoto 
et al. 1997) , and thus the luminosities in Table |^ are close 
to the total values. 

We plot the total X-ray luminosity versus Lb in Figure 
^ As expec ted for the relatively high Lx/Lb galaxies in 
our sample ( [Canizares et al. 1987 ), the slope is quite steep, 
d In Lx/d In Lb ~ 2, and is consistent with that predicted 
from steady-state cooling flow models (Nulsen et al. 1984; 



Sarazin 1997). Although our sample is too small to make a 
definitive statement, the slope appears to steepen for larger 
Lx and Lb. These galaxies include NGC 1399, NGC 5044, 
and NGC 5846 for which significant emission may be due to 
a diffuse intergroup medium. 

Also in Figure ^ we plot the emission from the soft and 
hard components separately. We take the emission of the 
soft components to be given by those of the cold component 
with temperature Tq. If Th < 2 keV then we add that emis- 
sion to the soft component as well. The hard components 
are those having Th ^ 5 keV. The slope of the soft emission 
appears to be slightly steeper than that of the total X-ray 
emission for the galaxies with smaller Lx. We can not dis- 
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0.001 



Figure 8. The ratio of (the luminosity of the cold compo- 
nent) to (the luminosity of the hot component) versus L^/L-g, 
where Lx = + ■ O^^ly galaxies with Tj] >~ 5 keV are plot- 
ted, while those with Tjj <~ 2 keV (indicating another hot gas 
phase) would appear somewhere in the upper right corner of the 
figure. All luminosities are evaluated in the 0.5-5 keV band. 



tinguish the slopes of the soft and hard components (though 
the galaxies with measured hard components do not extend 
to the highest Lx in our sample), but the hard component 
does appear to flatten somewhat at the largest Lb- 

In any event, the slope for the hard emission is still con- 
siderably steeper than 1 which indicates for the galaxies in 
our sample the hard components are not entirely associated 
with discrete sources. It is likely that with higher S/N data 
the temperatures of many of these hard components will be 
split into subcomponents some of which will have smaller 
temperatures consistent with emission from hot gas. 

As argued in previous studies of large samples of ellip- 
ticals with the Emstem Observatory (Canizares et al. 1987; 



Kim etl al. 1992) the ratio L^/ Lb is a good mdicator ot 
tke 



when the X-ray emission is dominated by hot gas (larger 
values) or by discrete sources (smaller values). In Figure ^ 
we plot the ratio of the luminosity of the cold component to 
the hot component, L'^ /L^ , versus L-^/Lb- Galaxies which 
have Th < 2 keV are not plotted since they have no mea- 
sured component presumably due to discrete sources, but 
would appear somewhere in the top right of the plot. We 
see an increase of L^/L^ with L^/Lb as expected, albeit 
with large scatter, and a slope slightly less than 1. Hence, 
from direct measurement of the ratio L J / L^ we have shown 
that as Z/x/Lb increases the proportion of the total emis- 
sion due to hot gas with respect to that of a hard compo- 
nent presumably due to discrete sources increases such that 
d\'a{L'^ / L^) /d\n{L^/ Lb) ~ 1. Interestingly, for a general 
slope a this relation implies that when L'^/Lb — > oo we 
have L^ — > Lb jo? which is consistent with expectations of 
a discrete component where L'^"'^ oc Lb- 




Figure 9. The temperature of the soft component versus the 
0.5-5 keV total X-ray luminosity. The soft temperature is equal 
to Tq for galaxies with Tfj >~ 5 keV and is the emission weighted 
value of Tq and Tjj for galaxies where Tjj <~ 2 keV. 



Finally, we plot the temperature of the soft component 
versus total X-ray luminosity in Figure ^ This temperature 
is taken to be Tq unless Th < 2 keV in which case it is 
set to the emission- weighted value of Tc and Th. Over the 
decade in luminosity Lx = (0.5 — 5) x 10*^ erg cm~'^ s~^ 
we see the expected rise of temperature of the hot gas with 
total X-ray luminosity. However, the trend disappears into 
scatter for larger Lx. For example, NGC 7619 and NGC 7626 
have low Tc for their X-ray luminosity. But they also have 
Th ~ 3 keV and thus their hot-gas temperatures may in fact 
be higher in reality. 

Another explanation is given by cooling flows. We infer 
mass deposition rates of Mgas > IAL0 yr"'^ for galaxies 
NGC 499, NGC 4636, NGC 7619, NGC 7626, and IC 4296 
which appear to have low soft temperatures f or their X-ray 



lumin osity. Hence, as in clusters of galaxies (Fabian et al 



1994) 



the presence of cooling flows may account for much 
of the scatter in the Lx — Tx relationship. 



4 THE EVIDENCE FOR EXCESS 

ABSORPTION AND ITS IMPLICATIONS 

For approximately half the galaxies in our sample the spec- 
tral fits are significantly improved when allowing for absorp- 
tion in excess of the Galactic value. This effect is most clearly 
illustrated for the galaxies where the cooling-fiow model pro- 
vides a substantially better fit than a single MEKAL (or 

In Figure ^ we display 



§3.1) 



Raymond-Smith) model (see 

the SISO spectrum of NGC 1399 and the best-fit cooling- 
fiow model obtained from jointly fitting the SISO and SISl 
data of the two available observations (see Table |^. When 
fixing A^'h to the Galactic value the best-fit model predicts 
more emission than observed for energies ~ 0.6 0.8 keV; 
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Figure 10. Fits of tlie absorbed cooling-flow model to NGC 1399 with fixed Galactic absorption (left) and with variable absorption 
(right). Only the SISO data is shown for sequence 80038000. 



= 575.0 for 257 dof giving Xicd = 2.24. These residuals 
around 0.8 keV cannot be eliminated by including another 
thermal component. However, if the absorbing column is al- 
lowed to be a free parameter then these residuals can be 
largely removed: = 341.3 for 256 dof giving Xrcd ~ 1-33. 
The resulting column density, A^h = l-8lo;2 ^ 10^^ cm~^, is 
over ten times the Galactic value of A^h = 0.13 x lO'^^ cm~^ 
( ^tark et al. 1992| ). 

Is this excess absorptio n an instrumental effect peculiar 
to ASCA? As explained in ^3.1, most previous fJOSAT stud- 



ies do not find evidence for excess absorption in early-type 
galaxies. This discrepancy between the spectral shape of 
ROSAT and ASCA below 1 keV has bee n observed, e.g. for 
the bright Seyfert gala xies MCG-6-30-15 (iNandra fc Pounds 
1992|)a|nd NGC 554 8 (llwasawa et al. 1997b| ) as well as the 
classic QSO 3C273 ( yaqoob et al. 1994). How ever, a recent 
analysis of MCG-6-30-15 by Orr et al. ( |l997a| ) finds that the 
~ 0.1 — 10 keV spectrum observed by Beppo-SAX (Parmai 
et al. 1<^9^ agrees extremely well with that of the ASCA SIS. 
Similar agreement between Beppo-SAX and ASCA has also 
been reported for 3C273 (Orr et al. 1997b). Furthermore, 



from a simultaneous ROSAT PSPC and ASCA observation 



of the Seyfert galaxy NGC 5548 ( [iwasawa et al. 1997b[ ) it 
has been found that the absorption determined from ASCA 
agrees with determinations from GINGA unlike the deter- 
mination from ROSAT. Given the agreement between the 
spectra of ASCA, SAX, and GINGA it would appear that 
the excess absorption measured by ASCA is not the result 
of an instrumental effect. 

We do note that the excess column density reported 
here, if due to a uniform screen of the galaxy, precludes 
any emission being detected in the carbon band of ROSAT. 
Of course, the excess absorption is likely to be distributed 
throughout the emitting gas, and the nearest unabsorbed 
parts of the galaxy will give detectable radiation in that 
band. (See Allen & Fabian (1997) for simulations of this 
effect for clusters of galaxies.) 



We mention that Rangarajan et al. (1995) did find ex- 
cess absorption for NGC 1399 by exploiting the spatial infor- 



mation provided by the ROSAT PSPC. By comparing the 
ratios of the spectrum of the inner region (r ^ 10 kpc) to a 
best-fit unabsorbed thermal model with the corresponding 
ratio of the spectrum of an outer annulus (r ~ 100 — 150 kpc) 
to its best-fit unabsorbed thermal model they clearly showed 
the diminution of photons with energies ^ 1 keV in the in- 
ner region (see their Figure 5). However, consistent with 
the general ASCA-ROSAT discrepancy stated above, excess 
absorption is not required when fitting absorbed (single - 
component) thermal models to the data (Jones et al. 1997). 

The excess A'h we have derived from the spectral fits is 
only the simplest result assuming a uniform absorber which 
reduces the X-ray emission according to an exponential form 
depending only on the energy. In this approximation the ex- 
cess columns A'^e ~ (1 — 3) x 10^^ cm~^ imply 

-,9 , 



atomic hydrogen, Mai 



5 X 10^ M, 



0' 



masses m 
when spread uni- 



formly over a 10 kpc radius. This level of absorber mass is 
very similar to the typical mass of the hot X-ray emitting 
gas in early- type galaxies, but Mabs is typically a factor of 
10 larger than the mass inferred from HI observations (e.g. 
Bregman, Roberts, & Giovanelli 1988; Bregman & Roberts 
1990). Similar disagreement is found for an absorpti on fea - 
ture in the field of NGC 4472 by Irwin & Sarazin ([l996|). 



However, these discrepancies with HI observations may be 
irrelevant since the absorber could be composed of molecu- 
lar gas (e.g. Ferland, Fabian, & Johnstone 1994). The mass 
of molecular gas can be measured with CO observations, but 
is typically factors of a few smaller than expected under the 
above simple assumptions (e.g. Irwin et al. 1997). However, 
excess absorbing mass above that predicted by the standard 
CO/H2 conversion has been reported from 158 micron ob - 
servations of the irregular galaxy IC 10 ( Madden et. al 1997 ). 
It should be emphasized that the uniform absorber model is 
certainly over-simplified and, e.g. a partial covering fraction 
and metallicity dependence could very well be required. 

The cold gas implied by the excess absorption is a nat- 
ural consequence of the multi-phase cooling-flow scenario 
(Nulsen 1986; Thomas et al. 1987; White & Sarazin 1987) 
wherein mass drops out of the flow at all radii where the 
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cooling time is less than the age of the galaxy (fa gc). I n par- 
ticular, for clusters of galaxies Allen & Fabian (L997) have 



demonstrated using color profiles of ROSAT PSPC data that 
clusters with tcooi ^ iagc have intrinsic absorbing material 
which is preferentially distributed in the core and of a quan- 
tity that is consistent with the amount of mass deposited by 
the cooling flow over tagc- 

Although the situation in early-type galaxies is likely 
complicated to some extent by heat sources, the success of 
the simple cooling-flow model over the single MEKAL model 
for the brightest galaxies and the associated required ex- 
cess absorption (see § [3.l| ) suggests that multi-phase cooling 
flows are operating in these galaxies. This is further sup- 
ported by the fact that these galaxies when fitted with two- 
temperatures indicate that they have a multi-phase medium 
(see §3.2); i.e. one component may represent the hot ambi- 
ent phase and the other the emission from the distributed 
mass drop-out component. In this scenario, the drop-out 
term should consist of a continuum of phases which will be- 
come evident when probed by detectors with better spectral 
resolution. 



5 CONCLUSIONS 

We present spectral analysis of ASCA data of 20 bright, 
mostly high L^/Lb, early-type galaxies (17 E, 3 SO) which 
is the largest sample of its kind reported to date. Fitting a 
single, absorbed, thin thermal model (MEKAL or Raymond- 
Smith) to the spectra fails to provide acceptable fits (e.g. 
Xred > 1-5) for 14 galaxies. Excess absorption significantly 
improves the fits for 10 of the galaxies. The abundances 
of the single-temperature models are very sub-solar, (Z) — 
0.19±0.12 Zq, in agreement with those reported by previous 
ROSAT and ASCA studies. 

Fitting instead an absorbed cooling-flow model yields 
fits of comparable quality to the MEKAL model for 11 galax- 
ies but is a significantly better fit for the other 9. The im- 
provement in the fits is dramatic for NGC 4472, NGC 5044, 
and IC 4296. The abundances obtained from the cooling- 
flow model, (Z) = 0.6 ± 0.5 Zq, significantly exceed those 
obtained with the MEKAL and Raymond-Smith models. 

We empirically tested the reliability of the MEKAL and 
Raymond-Smith plasma codes in the Fe-L band. In order to 
lessen the influence of other temperature components we re- 
stricted our analysis to 0.5-2 keV and performed fits over 
only the Fe-L band (taken to be 0.7-1.2 keV) and over the 
band excluding Fe-L (0.5-0.7 and 1.2-2 keV). We did not find 
large systematic trends expected if the plasma codes predict 
too much or too little Fe-L emission for a given tempera- 
ture and metallicity. However, unlike the MEKAL model, 
the Raymond-Smith model did give some evidence for these 
systematic trends which, as expected, indicates that the Fe- 
L region of Raymond-Smith is not as well modeled as with 
MEKAL. (The reasonable agreement of the abundances ob- 
tained by MEKAL and Raymond-Smith shows that such 
errors are not catastrophic.) 



Fit! 



Rayr 



Llug Lwo-Lenipt!i'aLure modelH (with MEKAL or 



-UmiLli mudt'ls) siguifiLauLl}' impiuved Lhe qualiLy 
of the fit in 16 out of 20 galaxies and gave formally accept- 
able values of x?cd ^ l-O- (We found that care needed to be 
exercised with XSPEC to achieve the global minima for 



the lowest S/N galaxies in the sample.) The galaxies having 
the highest S /N spectra in our sample (which coincident ally 
have the largest Lx/I/b) generally have a temperature for 
the cold component, Tq ~ 0.5 — 1 keV, and for the hot 
component, Th ~ 1 — 2 keV; i.e. both components are con- 
sistent with emission from hot ga f^. The lower S/N galax- 
ies (which coincident ally have the smallest L^/ Lb) typically 
have Tc 0.5 - 1 keV but Th ^ 5 keV. The temperature of 
this hot component is consistent with emission from discrete 
sources. The flux of this hot component is generally smaller 
than that of the cold component. 

In stark contrast to the single-temperature models, the 
abundances obtained from the two-temperature fits are ap- 
proximately solar; {Z) = 0.9 ± 0.7 Zq for two MEKAL 
models and (Z) = 0.7 ± 0.6 Zq for two Raymond-Smith 
models. The difference between the abundances determined 
from single-temperature and two-temperature models is not 
due to increased errors on the parameters but is an arti- 
fact of fitting intrinsically multi-temperature spectr a with 
single-temperature models (Buote & Canizares 1994). 

Our abundances disagree with the two-component mod- 
els fit to ASCA spectra of s everal of the galaxies in our 
sample by Matsumoto et al. (1997). We attribute the dis- 
crepancy to the different modeling procedures used and to 
the higher S /N data in our study owing to generally smaller 
extraction radii for the spectra. Possible abundance gradi- 
ents do not account for the discrepancy because we obtain 
approximately solar abundances even for some of the large 
S /N galaxies in our sample which have large extraction radii 
(e.g. NGC 1399, NGC 4472). The approximately solar abun- 
dances we find are consistent with multi-phase models for 
the evolution of the hot gas in ellipticals (e.g. Fujita et al. 
1997). However, though considerably larger than obtained 
in previous studies, these nearly solar abundances are still 
substantially less than the super-solar abundances predicted 
by the standard enrichment theories (e.g. Ciotti et al. 1991; 
also see Arimoto et al. 1997). 

Two-component models where one component is a 
cooling-flow model and the other a MEKAL give fits of com- 
parable quality to fits of two MEKAL models. Since the flux 
is shared by the two components, the inferred mass deposi- 
tion rates are smaller than when only a cooling-flow model 
is fit to the spectra. 

It should be emphasized that a radial temperature gra- 
dient in a single-phase gas unlikely accounts for the temper- 
atures seen in the large Lx/Lb galaxies. For example, NGC 



Recently Loewenstein & Mushotzky ( 1997 ) claim that ratios 



of the H-like to He-like Si lines in, e.g. NGC 1399 are only con- 
sistent with an isothermal gas with temperature near 1 keV and 



2. We 



inconsistent with the multi-phase solution we find in 
have examined such line ratios an d fi nd that the line ratios of 
the single-temperature models in S p.ll are very similar to those 
of the two-temperature models in j j3.2| (i.e. within ~ 15%). How- 
ever, caution must be exercised when computing these line ratios 
directly from the data because of the calibration errors due to 
the optical co nstants in the mirror around 2 keV (Gendreau & 
Yaqoob 1997| ). Since the key H-like Si line appears at 2.0 keV 



and for all of the galaxies in our sample this line is very noisy 
(e.g. Figure [lo[), and its intensity depends sensitively on the as- 
sumed continuum, we believe that any reliable constraints on the 
2.0 kcV H-likc Si line must await corrected response matrices. 
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4472 has roughly an equal contribution from phases with 
temperatures Tc = 0.76 keV and Th = 1.48 keV, but the 
temperature profile measured with the RO SAT PSPC varies 
fro m on ly ~ 0.75 — 1.1 keV over r < 4' (Irwin & Sarazin 
199( ^). ^ plausible physical explanation of the two phases is 
offered by a multi-phase cooling flow (Nulsen 1986; Thomas 
et al. 1987; White & Sarazin 1987) where one component 
represents the hot ambient phase and the other represents 
the emission from the distributed mass drop-out component. 
In this scenario, the drop-out term should consist of a con- 
tinuum of phases which should become evident when probed 
by detectors with better spectral resolution. 

We compared Lb to the luminosities of the hot and cold 
components determined by fits of two MEKAL models. Both 
the luminosities of the soft and hard components increase 
with Lb with a slope approximately 2, although there is 
some weak evidence that the cold component has a steeper 
slope than the hot component. For galaxies with Th J5; 5 keV 
we find that L J / ~ Lx /Lb ; i.e. there is a clear increase in 
the ratio of the luminosity of the cold component to the hot 
component as a function of Lx/Lb indicating that Lx/Lb 
is indeed a good indicator for the relative contribution of 
hot gas and discrete sources to the total X-ray emission. (A 
recent study of NGC 3923 by Buote & Canizares (1998) finds 
that the allowed contribution to the ROSAT X-ray emission 
from a population of discrete sources with Lx oc Lopticai is 
significantly less than indicated by L^ / measured in this 
paper if we associate oc Lopticai. This would support the 
conjecture that a significant fraction of the hard components 
measured for the low S/N galaxies is due to other phases in 
the hot gas.) 

Therefore, our investigation shows that the very sub- 
solar abundances obtained in previous studies are not the 
result of failures in the plasma codes (as proposed by Ari- 
moto et al. 1997). Rather, fitting multi-temperature models, 
without restricting the form of the models to a universal 
model for discrete sources, and setting extraction radii so 
that S/N is maximized, gives metal abundances near solar 
for both MEKAL and Raymond-Smith models. 

Moreover, fits to the X-ray spectra of half of the galaxies 
in our sample are significantly improved when allowing for 
absorption in excess of the Galactic value. Accurate deter- 
mination of the spatial distributions of the absorbing com- 
ponents and the various phases of the hot gas must await 
better quality data as must determination of the nature of 
the hard components found in the lower S/N galaxies in our 
sample. The next generation of X-ray satellites soon to be 
flown will have the capabilities to make significant progress 
in these areas. 
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